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^ ■ Abstract 

■ We consider the decoupling of neutrinos in the early Universe in presence of non- 

standard neutral current neutrino-electron interactions (NSI). We first discuss a 
Qh! semi-analytical approach to solve the relevant kinetic equations and then present the 

pi • results of fully numerical and momentum-dependent calculations, including flavor 

O , neutrino oscillations. We present our results in terms of both the effective number 

of neutrino species (A^eff ) a-nd the impact on the abundance of ^He produced during 
Big Bang Nucleosynthesis. We find that the presence of neutrino-electron NSI may 
^ ■ enhance the entropy transfer from electron-positron pairs into neutrinos instead of 

photons, up to a value of A^eff — 3.12 for NSI parameters within the ranges allowed 
by present laboratory data, which is almost three times the effect that appears for 
standard weak interactions. Thus non-standard neutrino-electron interactions do 
not essentially modify the density of relic neutrinos nor the bounds on neutrino 
properties from cosmological observables, such as their mass. 

Key words: Early Universe; Neutrinos; Non-equilibrium kinetics 



1 Introduction 



In the early Universe, neutrinos were kept in thermal contact with the elec- 
tromagnetic primordial plasma by rapid weak interactions with electrons and 
positrons. When the temperature dropped below a few MeV, these weak 
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processes became ineffective and the process of neutrino decoupling took 
place, while shortly after the pairs began to annihilate almost entirely 
into photons thus producing a difference between the temperatures of the 
relic photons and neutrinos. This difference can be easily calculated if we 
assume that neutrinos were completely decoupled when the pairs trans- 
ferred their entropy to photons, leading to the well-known temperature ratio 
Tr^/Ty = (11/4)^/'^ ~ 1.40102. Indeed, this simplified picture should be im- 
proved since some relic interactions between and neutrinos exist all along 
the annihilation stage, leading to a slightly smaller increase of the comov- 
ing photon temperature and to small distortions (at the percent level) of the 
neutrino momentum distributions. 

Presently, there exist compelling evidences for flavor neutrino oscillations from 
a variety of experimental data on solar, atmospheric, reactor and accelerator 
neutrinos (see e.g. [1,2]). These results are well understood by assuming that 
neutrinos have masses and mix, which in turn seems to point out the neces- 
sity of some new physics beyond the Standard Model (SM) of fundamental 
interactions. Interestingly, non-zero neutrino masses usually come with non- 
standard interactions (NSI) that might violate leptonic flavor and/or break 
weak universality. Recent analyses [3,4,5,6] have considered the neutral cur- 
rent NSI in a phenomenological way, showing that they can be bound using 
measurements of neutrino-electron scattering, as well as data from LEP and 
from related charged lepton processes. 

The aim of this paper is to study the neutrino decoupling process in presence of 
additional interactions between neutrinos and electrons, a possibility already 
noted in [3,5]. In this case, neutrinos could be kept in longer contact with 
and thus share a larger amount of the total entropy transfer than in the 
SM. Actually, if the non-standard neutrino-electron interactions were large 
enough, the neutrino momentum distribution would be significantly different 
from the standard case. In turn, this would modify the final yield of light nuclei 
during the epoch of Big Bang Nucleosynthesis (BBN), as well as the radiation 
content of the Universe, affecting the anisotropies of the Cosmic Microwave 
Background (CMB) and the power spectrum of Large Scale Structures (LSS). 
Our goal is to calculate how the decoupling is modified taking into account 
NSI with couplings which are still allowed by present laboratory data, and 
to discuss the possibility that cosmological observations can be used as a 
complementary way to bound these exotic scenarios. 

The paper is organized as follows. We begin in Sec. 2 by describing the for- 
malism adopted for the non-standard electron-neutrino interactions and sum- 
marize the current bounds from a variety of experimental data. We then con- 
sider the process of relic neutrino decoupling in the presence of non-standard 
electron-neutrino interactions, giving first an estimate by using a semi-analy- 
tical approach in Sec. 3. Finally, in Sec. 4 we report the results of the full 
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momentum-dependent numerical calculations for the neutrino spectra and 
the effect on the primordial ^He yield and other cosmological observables. 
We present our conclusions in Sec. 5. 



2 Non-standard neutrino-electron interactions 

The long-standing evidence of flavor change in atmospheric and solar neutrino 
experiments represents a strong indication of some new physics beyond the 
SM of fundamental interactions. For several years these experimental results 
have been interpreted in terms of neutrino masses (flavor oscillations) or by 
introducing new neutrino interactions. Recently, reactor and accelerator data 
have confirmed that neutrino masses are indeed the main mechanism explain- 
ing the atmospheric and solar anomalies, while the role of NSI can be only 
sub-leading. In any case, most extended particle physics models that account 
for neutrino masses naturally lead to new NSI, whose value strongly depends 
on the model. For instance, NSI may arise from the structure of the charged 
and neutral current weak interactions in seesaw- type extended models [7]. 

In the present analysis we will follow [3,5] and assume that new physics induces 
NSI only through the four fermion operators (z/z/)(//), where / is a charged 
lepton or quark, but not new charged lepton physics at tree level. In particular, 
since we are interested in the decoupling process of relic neutrinos, we consider 
only the NSI related to electrons which, together with the standard weak 
interactions, are described by the effective Lagrangian 

■^cff = -^^SM + X! '^NSI (1) 

which contains the four-fermion terms 



Cti = -2V2Gf E 5^/3 {uaYLu^Kn.Pe) (3) 
p 

for energies much smaller than the Z boson mass, as in our case, with Gp the 
Fermi constant and P = L, R = (l=P75)/2 the chiral projectors. Greek indices 
label lepton flavors {a, /5 = e, /i, r) and the Z couphngs are gi = — | + sin^ 9w 
and gR = sin^ 9w- 

The NSI parameters can induce a breaking of lepton universality {a = 
(3) or rather a flavor- changing contribution [a ^ (3). Their values can be 
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constrained by a variety of laboratory experiments, as discussed in [3,5,6]. In 
what follows we summarize the present bounds and refer the reader to the 
analyses mentioned above for all details. When considering the bounds on 
the parameters, it is important to notice that they are usually obtained 
taking only one-at-a-time, or at most combining two of them (such as the 
pair efe^) . This implies that the derived constraints are expected to be weaker 
but more robust if many NSI parameters are simultaneously included, since 
cancellations may occur. 

A remark is in order. Since isospin-changing weak interactions converting neu- 
trons into protons and vice versa are important for the BBN yields, one might 
wonder if it is legitimate to neglect the neutrino NSI with quarks. We argue 
that this approximation is well justified since, while relic neutrino distribu- 
tions affect n <-> p processes in a crucial way, vice versa is not true. The 
highly suppressed density of baryons with respect to electromagnetic parti- 
cles (rife/n^ ~ 6 X 10~^°) allows one to neglect completely neutrino scatter- 
ing on quarks for calculating the neutrino momentum spectrum. In addition, 
what contributes to the n ^ p rates are neutrino charged-current interactions, 
where possible non-standard terms are severely bounded by the accurate agree- 
ment of SM calculations for several processes with the data (e.g. leptonic and 
hadronic weak decays). For baryon thermalization, neutral current neutrino 
reactions are completely negligible with respect to electromagnetic interac- 
tions with and photons. Also note that any tiny non-standard effect in the 
weak rates is effectively taken into account by our prescription of rescaling the 
theoretically predicted neutron lifetime to the experimentally measured value 
(see [8] for details). 

Finally, we comment on a possible role of exotic four-fermion neutrino-neutrino 
interactions of the kind parameterized in Eq. (3), reporting the bounds quoted 
in [9]. If only left-handed neutrinos are involved, the accurate measurement 
of the Z— boson width constrains the NSI coupling to be at most of the same 
strength of the neutral current ones in the SM, and thus we expect that they 
will have a sub- leading effect on our results (see Sec. 4). Instead, for NSI terms 
coupling left-handed to right-handed components, a stringent BBN bound 
of O{10~^) applies, coming from the request that more than one neutrino- 
equivalent degree of freedom which was thermally populated in the early Uni- 
verse is excluded. Finally, extra interactions coupling only right-handed states 
might be large, but this case is of no interest for the process of neutrino de- 
couphng. 



2. 1 Bounds from tree level processes 
- Neutrino scattering experiments 
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The magnitude of the NSI parameters can be constrained from the analysis of 
data from neutrino-electron scattering experiments, which can probe the SM 
electroweak predictions with good precision (see e.g. [10]). 

First, we focus on z/g — e (data from the LSND experiment [11]) and z/g — e 
scattering (data from the Irvine [12] and MUNU [13] experiments). The total 
z/g — e cross section including NSI is 
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while the LSND measurement is 



:i.l7±0.17) ^ ' 



TT 



(5) 



In [5] this result was compared with the SM prediction a = l.OdGTGpmeE^/n — 
including electroweak corrections and best fit value for sin^ 9w — to obtain the 
following allowed regions around the SM value, assuming only one operator 
at a time 



- 0.07 (-0.14) < ef; < 0.11 (0.16) 

-0.99 (-1.14) < efg < 0.53(0.67) (6) 
\ei\<0A{0.5) 

kf.l<0.7(0.9) (7) 

at 90% CL (we added 99% CL bounds in parentheses). For more severe 
bounds are obtained using radiative effects, see later. In general, the allowed 
region in the e^,efg plane is an elliptic corona as given in Fig. 1, assuming 
^ea^ = for Q! 7^ e. As shown in [6], with the addition of data from i>g — e scat- 
tering the allowed LSND region can be substantially reduced. Since one must 
exchange Ql Qr with respect to Eq. (4) one obtains a perpendicular ellipse. 
The allowed regions for e^^^ are shown in Fig. 1, while the one-parameter 90% 
CL bounds are now [6] 



- 0.05 [-1.58] <ef'g < 0.12 [0.12] 
-0.04 [-0.61] < efg < 0.14 [0.15] 
<0.44 [0.85] 
< 0.27 [0.38] 
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Fig. 1. The ellipse corresponds to the values of (See^^ee) allowed at 99% CL by 
LSND data, while the dashed regions are those allowed at la and 99% CL when 
including data from antineutrino-electron experiments, as found in [6]. Here we 
assumed vanishing off-diagonal NSI parameters. We have calculated the full neutrino 
decoupling process for the indicated points (see Sec. 4). 

When leaving the corresponding flavor-changing NSI parameters free, the 
bounds are relaxed to the numbers in square parentheses: note the signifi- 
cant reduction especially for negative values. 

We also have data on z/^ — e scattering from the CHARM II collaboration [14]. 
The obtained results for the vector and axial vector e — couplings can be 
translated into the following 90% CL bounds on NSI parameters (see e.g. [5]) 



which imply that the y^ — e interactions must be very close to the SM predic- 
tions. 



As pointed out in [3], neutrino NSI can be also constrained by measuring the 



- 0.025 < e[„ < 0.03 
-0.027 < e± < 0.03 




(9) 
(10) 



LEP data on e~^e 



z/z/7 
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e~^e~ — * z/z/7 cross section. This is actually the only way to get bounds on 
the e:^;.^ parameters from laboratory data. The approximate limits for these 
parameters can be extracted from Fig. 3 of [4], 



0.7 < < 0.5 



-0.5<5«<0.6 (11) 

at 99% CL when the other parameter is left free. For the z/g e NSI the bounds 
are comparable to those obtained from neutrino-electron scattering. 



2.2 One loop effects 



The non-standard neutrino-electron interactions also lead to corrections at the 
one-loop level to processes such as the decays of the electroweak gauge bosons 
or lepton flavor violating decays of charged leptons. These corrections arise 
from effective non-renormalizable interactions and in principle their computa- 
tion requires the knowledge of the complete theory leading to these effective 
terms in the low-energy regime. Nevertheless for A 3> mw, where A is the 
energy scale setting the limit of validity of the effective theory, the leading 
term is independent of the specific theory and has a logarithmic behavior, 
ln(A/miy) [5]. Using conservatively ln{A/m\y) ~ 1, one gets the following 
90% CL bounds from the decay rates of the electroweak gauge bosons 

\e^f\ < 0.5 (12) 



which are of the same order of Eq. (11). 

Similarly, from the strong experimental limit on the branching ratio Br{fi~ 
e~e~^e~) < 10~^^, it is possible to obtain a severe bound on the flavor changing 
parameter (90% CL) 

kfJ<5xlO-^ (13) 



2.3 Summary of bounds on NSI parameters 



After reviewing the bounds on NSI parameters, we conclude that they 

• restrict the £^ as shown in Fig. 1; 

• force the e^^ to be very small and negligible for the purposes of the present 
paper; 

• allow large e^^, almost of order unity, see (11); 
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• as far as flavor changing parameters are concerned, allow large e^, (7), 
moderate e^^, (10), and very tiny e^, the latter being too small to have any 
impact on neutrino decoupling dynamics. 

Again it is useful to remind the reader that the bounds on NSI parameters 
are relaxed when more than one parameter is allowed to vary simultaneously, 
as shown in [6]. 

As noted in [5] , the leptonic measurement of sin^ 9\y at a neutrino factory will 
be sensitive to values of the i^e,/i ~ e NSI parameters of O{10~^). Additional 
information on the NSI can be extracted from a consistency comparison of the 
results of solar neutrino experiments with those of an experiment detecting re- 
actor neutrinos such as KamLAND. The former depend on the NSI parameters 
through the matter potential felt by solar neutrinos and/or via the neutrino 
neutral current detection. Instead, KamLAND data is essentially unaffected 
by the NSI. Ref. [5] found that the combination of SNO and KamLAND data ^ 
will provide bounds on the e^^ comparable but slightly better to those in Eqs. 
(11) or (12). However, it has been noted in several works [15,16,17,18,19] that 
while NSI are expected to play a subdominant role on the oscillations of atmo- 
spheric and solar neutrinos, there exist degenerate directions in the parameter 
space where large NSI are still allowed. Future data from short and medium 
baseline neutrino beams will help to resolve these degeneracies [20]. 



3 Neutrino decoupling in presence of NSI 

3.1 Delayed neutrino decoupling in the instantaneous limit 

In the early Universe, thermally produced neutrinos were in equilibrium with 
other particles down to temperatures of few MeV, when weak interactions 
became ineffective and neutrino decoupled from the plasma. As a first ap- 
proximation, the thermal Fermi-Dirac momentum spectrum is preserved after 
the freeze-out of weak interactions, since neutrinos decoupled when ultra- 
relativistic and both neutrino momenta and temperature redshift identically 
with the Universe expansion. A neutrino chemical potential /ij, would exist in 
the presence of a neutrino- antineutrino asymmetry, but it was shown in [21] 
that the stringent BBN bounds on /ij,^ (for an updated limit see e.g. [22]) ap- 
ply to all flavors, since neutrino oscillations lead to flavor equilibrium before 
BBN. Thus we shall ignore a relic neutrino asymmetry in the following. 

It proves useful to define the following dimensionless variables instead of time, 

^ The Borexino detector was considered in [4]. 
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momenta and photon temperature 

x = mR y = pR z = T^R, (14) 



where m is an arbitrary mass scale which we choose as the electron mass rrie 
and R is the Universe scale factor. The function R is normalized without loss 
of generality so that R{t) 1/T at large temperatures, T being the common 
temperature of the particles in equilibrium far from any entropy-transfer pro- 
cess. With this choice, R~^ can be identified with the temperature of neutrinos 
in the limit of instantaneous decoupling. After decoupling, neutrinos enter cos- 
mological observables mainly via their energy density. Since they also remain 
relativistic for most of the cosmological evolution — down to a temperature of 
the order of their sub-eV mass scale — it is customary to parameterize their 
contribution to the radiation energy density Pr in terms of the effective num- 
ber of neutrinos N^ff [23,24] 



where Zq{x) describes the photon to neutrino temperature ratio in the limit 
of instantaneous neutrino decoupling. For x ^ 1 the annihilation phase is 
over, and zq — *• (11/4)^/^ ~ 1.4010. The energy densities p° and refer re- 
spectively to the photon plasma and to a single neutrino species in the limit of 
instantaneous decoupling, while p^ is the actual energy content of the photon 
plasma and Pu+x the total energy content of weakly interacting particles (in- 
cluding possible exotic contributions). The second equality in Eq. (15) follows 
when only the three active neutrinos contribute to Pu+x', eventually, the actual 
photon temperature evolution accounts for the second factor in Eq. (15) and 
the possible energy- density distortion in the a-th neutrino flavor is given by 
6a = {pu^ - pD/pI- Note that from Eq. (15) it follows 

PR = (^1 + ^4 ^efr) P7 , (16) 



which, replacing zq with (11/4)^/^, is often used in the literature to define 
A^cff in the asymptotic limit x ^ 1. Clearly, well after e"^ annihilation, three 
thermally distributed neutrinos correspond to A^cfr = 3 in the instantaneous 
decoupling limit. 

In order to estimate analytically how large would be the impact of NSI on 
the decoupling mechanism, let us discuss a simple toy-model, which we shall 
later compare with the results obtained solving the relevant kinetic equations. 
Consider the general case where 3 — A^' neutrinos interact via standard weak 
processes, while the remaining A^' have extra contributions from NSI which 
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enhance the interaction rates over e^. In the instantaneous decouphng hmit for 
both the neutrino species, the plasma can be described by a single additional 
parameter, the non-standard neutrino temperature T' or, equivalently, by w = 
T' R>1. From Eq. (15), one obtains 



eflf 



{3-N')+N'w' 



3 + N\w' 



(17) 



where the second step reflects the fact that, for a single NSI neutrino species, 
the distortion is given by 6' = {w^ — 1)- If we denote with Xd = TdR the epoch 
of the standard neutrino decoupling, and parameterize with x'^ = T^i? the 
decoupling temperature of the non-standard species, we can easily calculate 
the functions w{x) and z{x) from the conservation of entropy per comoving 
volume, 



'W[X) 

w{x) 
w{x) 
z{x) 



z{x) = 
z{x) = 

-- Wd = 



N's^{xj^)+Sj(xa)+s^± (xg) 
N's^(x'^)+s-,{x'^)+s^±{x'^) 



Wd 



S7(2^d) + «e±(^d) 



s^{x)+s^±{x) 



1/3 



1/3 



1/3 ^ 



X < Xd, 
Xd<X< x'^, 



X > x' 



di 



where s^,, and Se± are respectively the speciflc entropy of one u-u species, 
of the photons and of the e^. 

In Fig. 2 we show the results for N^.^ in the instantaneous decoupling limit 
as function of the decoupling temperature and for various choices of A^'. 
As expected, iVeff grows if neutrinos decouple at a smaller temperature and 
shows two regimes. For very small NSI contribution to the neutrino interaction 
rates we have Wd — ^ and N^^i ^ 3 since neutrinos decouple long before any 
annihilation. On the other hand, for A^' neutrinos tightly coupled to the 
electromagnetic component down to very small temperatures when all pairs 
have already annihilated, 

/22 + 7A^'\^^^ , , 

z ~ — -— , ii? ~ 2;. (18) 



Thus, in the tightly coupled limit, one flnds N^s — > {7.0,9.7, 11.5} for A^' = 
{1, 2, 3} respectively, as shown in Fig. 2. Note that one obtains a contribution 
to the radiation energy density of order N^s ~ 4 if neutrinos are kept in 
equilibrium down to temperatures of 0.2 — 0.3 MeV. 
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Fig. 2. Contribution to the radiation energy density parameterized in terms of N^s, 
when 1, 2 or 3 neutrino species are coupled to the electromagnetic plasma until a 
given value of the decoupling temperature T^. 



3.2 Decoupling temperature with NSI 



In order to dynamically predict the decoupling temperature T'^ when NSI are 
present, we have to solve the relevant kinetic equations taking into account the 
new couplings e. For a first analytical estimate, we can compare the interaction 
rate F with with the expansion rate of the Universe given by the Hubble 
parameter H. Since H oc T^, and F ~ (1 +e)^T^, the decoupling temperature 
at which H = T decreases as ~ (1 + . As we have argued in the 

previous section, in order to produce changes of 0{1) in A'efr the decoupling 
temperature should be lowered down to ~0.2-0.3 MeV, i.e. should be one 
order of magnitude smaller than for ordinary neutrinos. This implies e > 20, 
which would largely exceed present laboratory bounds. We thus expect at most 
changes in A'^efr of (9(0.01-0.1). To improve further our treatment, we can use 
an approximate solution of the kinetic equation for the neutrino distribution 
function /^(x, y), 



Hx 



dx 



rcoU 



(19) 
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Fig. 3. Contours of equal neutrino decoupling temperature in MeV for different 
values of the parameters characterizing the Ve — ^ and the Vr — ^ diagonal NSI. 

where /^°" is the collision term that includes all relevant neutrino interaction 
processes. Following [25], we shall keep only the direct reaction term in the 
collision integral using the appropriate matrix elements for neutrino elastic 
scattering and the inverse annihilation e^e~ vv in the relativistic limit 
(me — »• 0). Assuming that all particles that interact are close to equilibrium, 
the collision integral can be estimated in the Boltzmann approximation leading 
to 



with g\ — > (1 + for the case of z/g. This equation can be easily integrated 
in X, leading to the values T^iye) — 1-8 MeV and Td{i'^^r) — 3.1 MeV for 
an average momentum of (y) = 3.15. This exercise can be repeated in the 
presence of neutrino-electron NSI in order to find their influence over the de- 
coupling of neutrinos, adding e^^^ or e^^^- to the couplings and gji. The 
calculated decoupling temperature is shown in Fig. 3 for the case of non-zero 
^g^^ and e!^^^. In general, a significant increase of the NSI parameters from 
the SM prediction leads to a larger interaction of neutrinos with and thus 
to a lower decoupling temperature. However, for a region close to the pair of 
values {e^,e^) that minimize the interaction by accidental cancellation with 
SM couplings, the decoupling temperature is significantly raised. Finally, one 
can relate the value of the decoupling temperature in presence of NSI with 
an estimate of the change in A^efr following the instantaneous decoupling ap- 
proximation previously described. The results are shown in Fig. 4. As already 
noticed, large modifications of A^efr require that the neutrino-electron interac- 
tions are much larger than ordinary weak processes and are thus excluded by 
laboratory bounds. 
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Fig. 4. Contours of equal Ncg for different values of the NSI parameters as found in 
the instantaneous decoupling approximation. 

Nonetheless, this study suggests that additional distortions in the neutrino 
spectra comparable to or larger than the ones predicted in the SM are possible, 
and it is still interesting to assess exactly how large these effects might be. 
However, the treatment followed till now assumes a series of approximations 
(such as rrie — >^ 0, no QED corrections to the plasma, no effect of the active 
neutrino oscillations, etc.) that in the next section we shall check using a fully 
numerical and momentum-dependent calculation. 



4 Numerical calculation: momentum-dependent kinetic equations 



4.1 Method 



The instantaneous decoupling approximation does not take into account that 
neutrino-electron interactions are more efficient for larger neutrino energies. 
Thus one expects a priori that more energetic neutrinos will remain longer in 
thermal contact, leading in general to non-thermal distortions in the neutrino 
spectra and a slightly smaller increase of the comoving photon temperature z. 
This is actually the case for standard weak interactions, as noted in previous 
works (for early references, see [26] and the full list given in the review [25]). 

A proper calculation of the process of non-instantaneous neutrino decoupling 
demands solving the momentum-dependent Boltzmann equations for the neu- 
trino spectra, a set of integro-differential kinetic equations that are difficult 
to solve numerically. In the early 1990s several works [27,28,29] performed 
momentum-dependent calculations assuming some approximations, such as 
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Boltzmann statistics for neutrinos, while the full numerical computation was 
later carried out in [30,31,32,33]. Finally, a further refinement involves the in- 
clusion of finite temperature QED corrections to the electromagnetic plasma, 
[34,35], while the last works also include the effect of flavor neutrino oscilla- 
tions [36,37]. 

Here we will follow our previous work [37]. In particular, we will describe the 
neutrino ensemble in the usual way by generalized occupation numbers, i.e. by 
3x3 density matrices q for neutrinos and anti- neutrinos as described in [38,39], 
with elements Qap where = e,yU, r. The diagonal elements correspond to 
the usual occupation numbers of the different flavors, while the off-diagonal 
terms are non-zero in the presence of neutrino mixing or flavor-changing NSI. 
The equations of motion for the density matrices relevant for our situation of 
interest in an expanding Universe are 



i {dt - Hp dp) Qp 



8V2Gfp\ 



C[0p\ , (21) 



where mw is the W boson mass. We use the notation Qp = g{p,t) and [■, ■] 
denotes the commutator. The vacuum oscillation term is proportional to M^, 
the mass-squared matrix in the flavor basis that is related to the diagonal one 
in the mass basis diag(m^, mg, m|) via the neutrino mixing matrix U, which 
in turn depends on the neutrino mixing angles 9i2, 023 and ^13 (we assume CP 
conservation). From a global analysis of experimental data on flavor neutrino 
oscillations, the values of mixing parameters can be extracted. As a reference, 
we take the best-fit values from [1] 

10-5 eV^ ' 10-3 eV^ '^'^'^^3.sL 1 = (8.1,2.2,0.3,0.5,0) (22) 



where s^^- = sin^ 9ij for i,j = 1,2, 3. In Ref. [37] we have shown that the results 
are not essentially modified for other values within the allowed regions, in 
particular for non-zero ^13 close to the experimental upper bound. Since we 
assume maximal 623 and zero 613 the spectra of z/^ and will be the same 
after decoupling. 

The effects of the medium in Eq. (21) are contained in the collision term C[-] 
and in a refractive term that corresponds in Eq. (21) to the term proportional 
to the diagonal matrix E, that represents the energy densities of charged 
leptons. See [37] for a description of both terms. The kinetic equations (21) 
should be numerically solved along with the equation governing the evolution 
of the electromagnetic temperature during the process of annihilations. 
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This can be found from the continuity equation for the total energy density p, 
^ = -3H (p + P), (23) 



with P the pressure, which can be cast into a first order differential equation 
for z as function of x, see e.g. 



We have generalized the results of [37] calculating the full evolution of neutrino 
during decoupling in presence of non-standard interactions with electrons, 
which modify the equations in two different ways. First, the couplings are 
modified inserting the parameters e^^^, e^'^ or , which implies a variation 
on the coUisional term in Eq. (21). Second, neutrino refraction in the medium 
is modified by the NSI, which we took into account in the refractive term in 
Eq. (21) 




i^matt = -T^P- (24) 



where pe is the energy density of the background electrons and positrons and 
£af3 = ^a/3+^f/3! siucc matter effects are sensitive only to the vector component 
of the interaction. The above equation includes only those NSI parameters that 
are not constrained to be smaller than (9(0.1) (see Sec. 2), and thus represents 
the leading-order correction to the refractive index. 

We have numerically solved the kinetic equations Eq. (21) using a discretiza- 
tion in a grid of 100 dimensionless neutrino momenta in the range yi G 
[0.02,20]. We start to compute the evolution of the system at a value Xin, 
when weak interactions were effective enough to keep neutrinos in equilib- 
rium with the electromagnetic plasma but flavor oscillations are suppressed 
by medium effects. The system of equations is then solved from Xin until a 
value of X when both the neutrino distortions and the comoving photon tem- 
perature z are frozen, approximately at Xfin — 35. In some cases, we made the 
calculation for different initial conditions and concluded that the results are 
stable within per cent accuracy unless we start the numerical evaluation at 
Xin ^ 0.2. For other technical details, we refer the reader to [37]. 



4-2 Numerical results 



We have numerically calculated the evolution of the neutrino density matrix 
solving the system of Eqs. (21) and (23), during the full process of neutrino 
decoupling. Since a complete scan of all possible combinations of the NSI 
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0.1 



1 

X = m„ R 



10 



Fig. 5. Evolution of the distortion of the t'e and v^^r spectra for a particular co- 
moving momentum (y = 10) with standard weak interactions (solid line) and with 
Ur — e NSI (dotted line). The line labelled with Ty corresponds to the distribution 
of a neutrino in full thermal contact with the electromagnetic plasma. 

parameters can not be done, we have performed the calculation for a selection 
of values for the NSI parameters, representative of the regions allowed by 
present bounds (see Sec. 2). 

As an example, in Fig. 5 we plot the distortion of the neutrino distribution as a 
function of x for a particular neutrino momentum [y = 10), both in the case of 
SM weak interactions and in presence of large NSI between tau neutrinos and 
e^, corresponding to e!^^ = e^^ = 2. The behavior of fi, has been previously 
described e.g. in [31,33,37]. At large temperatures or x < 0.2, neutrinos are 
in good thermal contact with and their distributions only change keeping 
an equilibrium shape with the photon temperature [exp{y/z{x)) + 1]~^ (the 
line in the figure). In the intermediate region 0.2 < x < 4, the standard 
weak interactions become less effective in a momentum-dependent way, lead- 
ing to distortions in the neutrino spectra which are larger for i/g's than for 
the other flavors. Finally, at larger values of x neutrino decoupling is complete 
and the distortions reach their asymptotic values. In the case of large u^- — e 
NSI, neutrinos are kept in thermal contact with for a longer time, leading 
to larger distortions for all neutrino flavors, in particular for u^ .,- (although 
Sg^-^ = e^'^ = 0, the electron and muon neutrinos also acquire significant 
distortions due to the effect of flavor oscillations). For the particular neutrino 
momentum in Fig. 5, if e!^^ = e^^ = 2 the final values of the distribution are 
9.8% for the v^s and 14.3% for the z/^,r's larger than in the limit of neutrino 
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Table 1 

Frozen values of Zf^^, the neutrino energy densities distortion 5a, N^s and AYp in 
presence of f e — e NSI. 





r-R 


2fin 


0e{70) 




AT 


AYp 




U 





l.o97o 


0.73 


0.52 


3.046 


2.1x 10 


-4 


0.1 


0.5 


1.3969 


1.20 


0.72 


3.062 


2.8x10" 


-4 


-1.0 


1.0 


1.3966 


1.32 


0.77 


3.067 


3.0x10-^ 


-4 


-1.5 


0.1 


1.3976 


0.81 


0.56 


3.048 


2.0x10" 


-4 


-1.5 


-1.0 


1.3969 


1.17 


0.71 


3.061 


2.7x10" 


-4 


-1.0 


-1.5 


1.3977 


1.45 


0.83 


3.060 


2.9x10" 


-4 


-0.07 


-1.0 


1.3973 


0.99 


0.64 


3.055 


2.6x10" 


-4 


4.0 


4.0 


1.3812 


9.47 


3.83 


3.357 


11.7x10 


-4 



decoupling before any e annihilation. For comparison, the corresponding val- 
ues for standard weak interactions are 4.4% for the u^s and 2% for the v^^tS. 

As in Ref. [37], we will summarize the results in terms of the following frozen 
values: dimensionless photon temperature Zfin, the fractional changes in the 
neutrino energy densities 5a and the asymptotic effective number of neutri- 
nos A'^eflf as defined in Eq. (15). We have considered four main sets of NSI 
parameters: 

(1) Only ^ee^- results are summarized in Table 1. We have calculated var- 
ious combinations of the two Ve ~ & NSI parameters, that correspond to the 
points shown in Fig. 1. For illustrative purposes, we considered the case with 
extremely large NSI = = 4 (already excluded by laboratory bounds), 
where the distortions in the distribution functions of neutrinos lead to a sig- 
nificant change in Nf,^. 

(2) Only e^:^: our results are summarized in Table 2 for values of both pa- 
rameters up to order unity in absolute value. Again, we also show one case 
with larger NSI e^^ = e^^ = 2, disfavored by terrestrial experiments. 

(3) Combinations of e^^^^ and e:^;.^: our results are summarized in Table 3. Here 
we varied the values of the e^^^ pair, while the — e NSI parameters were 
fixed to a value close to the current upper bound from laboratory experiments, 
e^, = e^^ = 0.5. 

(4) Cases where the flavor- changing NSI parameters e^^^^ are non-zero, that for 
simplicity we consider real. We consider both cases with vanishing diagonal 
NSI and with all possible NSI parameters with the largest allowed values 
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Table 2 

Same as Table 1 but for — e NSI. 









6e{%) 




iVcff 












1.3978 


0.73 


0.52 


3.046 


2.1x10" 


-4 


0.5 


0.5 


1.3970 


0.90 


0.79 


3.059 


4.0x10" 


-4 


1.0 


1.0 


1.3966 


1.24 


1.43 


3.079 


4.8x10" 


-4 


-1.0 


-1.0 


1.3952 


1.24 


1.45 


3.092 


5.3x10" 


-4 


2.0 


2.0 


1.3911 


2.04 


2.96 


3.168 


10.0x10 


-4 



Table 3 

Same as Table 1 when both the v^. — e and Vr — ^ NSI are non-zero. In all the cases 
we have fixed el^^ = = 0.5. 







Zfin 












0.1 


0.5 


1.3963 


1.30 


0.95 


3.073 


3.6x10" 


-4 


-1.0 


1.0 


1.3960 


1.41 


0.99 


3.077 


3.7x10" 


-4 


-1.5 


-1.0 


1.3963 


1.28 


0.94 


3.073 


3.4x10" 


-4 


-1.0 





1.3977 


0.54 


0.66 


3.047 


2.6x10" 


-4 



Table 4 

Same as Table 1 when the flavor-changing NSI parameters Ser^ non-zero. The 
last two cases correspond to the results for all NSI parameters non-zero with values 
given by the bounds in Eqs. (8). 















2fin 


5e{%) 






AYp 
















0.4 


0.7 


1.3960 


1.32 


1.03 


3.077 


3.9x10" 


-4 














0.85 


0.38 


1.3956 


1.46 


1.16 


3.085 


4.4x10" 


-4 


0.12 


0.15 


0.5 


0.5 


0.85 


0.38 


1.3949 


1.68 


1.38 


3.098 


5.0x10" 


-4 


-0.61 


-1.58 


-0.5 


-0.5 


-0.85 


-0.38 


1.3948 


1.69 


1.41 


3.100 


5.1x10" 


-4 


0.12 


-1.58 


-0.5 


0.5 


-0.85 


0.38 


1.3937 


2.21 


1.66 


3.120 


6.0x10" 


-4 



according to Eqs. (8). 

The last column [AYp) in Tables 1-4 reports the variation in the mass yield 
of ^He synthesized during BBN, which predicts typically Yp ~ 0.248 for the 
baryon density = 0.0223 ±0.0008 favoured by CMB anisotropy studies [40]. 
Neutrinos affect the outcome of BBN in a flavor independent way through their 
contribution to the radiation energy density (change in Ngs). In addition, the 
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role played by electron neutrinos and antineutrinos in the weak reactions that 
convert neutrons and protons is very important for fixing the primordial pro- 
duction of '^He, which is thus the nucleus mostly sensitive to non-standard 
physics in the neutrino sector. To predict quantitatively the change Al^ in- 
duced by NSI, we take exactly into account the modified thermodinamical 
quantities, while treating perturbatively the change in the weak rates. This is 
justified given the relatively small effects we are considering. We closely follow 
the treatment given in [37], to which we refer the reader for further details. 



The results are generally in agreement with the approximate ones predicted 
in Sec. 3 (see in particular Fig. 4), apart for minor differences due to the 
improved treatment. Also, as one naively expects, large e^^^^ enhance mostly 
the distortions in the electronic flavor, while large e!;^^^ lead to larger distortions 
in the tau neutrino spectra. However, the differences between the spectra of 
the three neutrino flavors are reduced by the partial re-shuffling of the entropy 
transfer due to the effects of neutrino oscillations. 



In general, we flnd that the presence of non-zero NSI enhances the transfer of 
entropy from to neutrinos, leading to values up to A^^efr — 3.12 when all NSI 
parameters are close to the boundaries of the allowed regions discussed in Sec. 
2, which is almost three times the departure of A'efr from 3 that already exists 
for SM weak interactions (the flrst row in Table 1). We also obtained some 
signiflcant departures of iVeff from the standard value, such as iVeff ~ 3.36 for 
= = 4, but only for NSI parameters well beyond the allowed regions. 
On the other hand, for particular choices of the diagonal NSI parameters 
that minimize the u — e interaction we flnd a reduction in the flnal value of 
the distortions or N^s with respect to the SM value, as expected from our 
discussion in Sec. 3.2. 



Finally, we mentioned in Sec. 2 that neutrino-neutrino interactions are poorly 
constrained. In the case of interactions involving left-handed neutrinos only, 
they can still be twice as effective as in the SM, as reported in [9]. Since en- 
hanced neutrino-neutrino interactions will redistribute the distortions more ef- 
flciently, they can modify the results described in this section. In order to check 
this, we doubled the intensity of the contributions of the neutrino-neutrino in- 
teractions to the collision integral and repeated the calculations for some of 
the cases considered above. We found that the results are basically unchanged, 
conflrming the unimportant role played by neutrino-neutrino collisions in the 
process of neutrino decoupling. 
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5 Conclusions 



The process of relic neutrino decoupling in the early Universe is sensitive to 
the strength of the interactions between neutrinos and the plasma formed 
by electrons and positrons. If neutrinos were kept in longer thermal contact 
with them than in the standard case, they would share a larger fraction of 
the entropy release from annihilations. This would affect the predicted 
characteristics of the cosmic background of neutrinos (CuB), which in turn 
could modify the late evolution of the Universe and the bounds on neutrino 
properties from the analysis of cosmological observables. 

In this paper we have considered how the decoupling of relic neutrinos is modi- 
fied in presence of non-standard neutral current neutrino-electron interactions. 
First, we have provided a rough estimate of the size of NSl couplings needed to 
affect in a relevant way the properties of the Cz/B and of related observables. 
We find that needed couplings are by far larger than the existing laboratory 
bounds, as also confirmed by a semi-analytical solution of the relevant kinetic 
equations. Thus, NSl can only play a minor role in shaping the Ci^B. However, 
NSI might still contribute with comparable or larger distortions in the neutrino 
spectra than predicted in the SM. In order to quantify these deviations, we 
have performed fully numerical and momentum- dependent calculations of the 
density-matrix equations relevant for neutrino evolution in the early Universe, 
including the effects of flavor neutrino oscillations. 

Typically, we find that the enhancement in observables like iVeff, which mea- 
sures the change in the radiation content, or the ^He mass yield in the BEN 
can be up to three times larger than the ones found in the standard case 
[37], for values of all NSl parameters close to the limits placed by laboratory 
experiments. Nevertheless, even the variation of up to 0.2% that we found 
in the ^He abundance is yet too small (about one order of magnitude) when 
compared to the observational error on Yp, which is unfortunately dominated 
by systematics, and as such difficult to pin down in the near future. Instead, a 
value of A'^efj ~ 3.12, which is still possible within the present parameter space 
for NSl couplings, might even be barely detectable. Indeed, it has been shown 
that the CMB satellite PLANCK will soon provide temperature and polariza- 
tion data that will probably measure the value of N^s with an uncertainty of 
order a{Naa) — 0.2 [41,42] (for a previous, more optimistic forecast see [43]), 
or better when combined with data from a large galaxy redshift survey such 
as SDSS [44]. Future CMB missions may reach the sensitivity of 0.04 — 0.05 
needed to test the standard scenario [42], and might achieve a 2 a hint for 
Ncs — 3.12 and eventually test the effect of large neutrino-electron NSl. 

However, barring some extreme cases, we conclude that the prediction of 
A'^eff — 3 within a few % is quite robust even when taking into account neutrino- 
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electron NSI of the four-fermion type. Thus, their existence can not modify 
in a significant way the bounds on neutrino properties from cosmological ob- 
servables, in particular on their masses, as recently reviewed in [45]. Since it 
is likely that future experiments may narrow further the allowed range of NSI 
couplings, their role in a cosmological context will be even smaller, at most 
a sub-leading correction to the standard prediction of the Cz/B properties. 
Turning the argument around, however, one can conclude that a significant 
deviation of N^s from 3 may require major revisions of the cosmological model, 
like the introduction of new relativistic relics and/or of an exotic thermal his- 
tory, or both (see e.g. [46,47]). 
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